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Abstract 

We extend our new approach for numeric evaluation of Feynman diagrams to integrals that include fermionic 
and vector propagators. In this initial discussion we begin by deriving the Sine function representation for the 
propagators of spin-i and spin-1 fields and exploring their properties. We show that the attributes of the spin-0 
propagator which allowed us to derive the Sine function representation for scalar field Feynman integrals are shared 
by fields with non-zero spin. We then investigate the application of the Sine function representation to simple QED 
diagrams, including first order corrections to the propagators and the vertex. 
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1 Introduction 

In conjunction with several colleagues, four of the present 
authors are in the process of developing a novel framework 
for non-perturbative numerical calculations in Quantum 
Field Theory. This Source Galerkin Method g § |, § 
is free of many of the problems that plague Monte Carlo 
calculations, in part because it works on the continuum 
and treats bosons and fermions in a largely symmetrical 
manner. A useful and surprising by-product of this "ex- 
act" method is that some of the computational techniques 
we developed in order to implement it are also applica- 
ble to the numerical evaluation of Feynman diagrams in 
conventional perturbation theory. 

The resulting new method of numerically evaluating 
Feynman diagrams, called the Sine function representa- 
tion, was introduced by three of us in || , and extensively 
tested in [Q . It relies on the observation that almost any 
integral can be approximated as an infinite sum, via the 
use of generalized Sine functions.^] Using this representa- 
tion, we develop approximate versions of the field propa- 
gators where the spatial (or momentum) dependence only 
appears in terms like exp(— x 2 ), so that all the integrals 
over vertex locations or internal momenta are reduced to 
Gaussian integrals, which can be performed analytically. 
After these integrals are evaluated, the diagrams are rep- 



1 We use the capital S in Sine to differentiate these functions from 
the standard form, sinc(x) = s'mx/x. 



resented as a multi-dimensional but rapidly convergent 
sum, which is evaluated numerically. 

The initial investigation of the Sine function represen- 
tation of the propagators of fields with non-zero spin 4] 
and their application to Feynman diagram calculations | £ ] 
strongly suggested that the Sine function representation 
would be useful for the calculation of diagrams in gauge 
theories. However, until now our systematic development 
of the Sine function representation has been limited to 
scalar fields, which allowed us to develop our procedures 
in a relatively simple, index-free environment. Here we 
demonstrate that our evaluation methods generalize di- 
rectly to fields with spin, and thus provide a powerful 
tool for calculations in more general theories. In this pa- 
per we focus on diagrams whose analytic properties are 
well understood, and calculate a variety of second order 
QED diagrams. We will extend our calculations to more 
complicated graphs in future work. We argue that these 
calculations — in conjunction with the higher order evalu- 
ations we have performed for scalar diagrams — provide 
strong evidence that the Sine function representation will 
be a useful calculational tool for fermionic and gauge the- 
ories. 

In the next section we recapitulate the most impor- 
tant facts about the generalized Sine function and in the 
following section we derive the approximate versions of 
vector and Dirac field propagators, checking their accu- 
racy by comparing them to the exact forms. In Section 
4 we explicitly evaluate several diagrams, showing that 
we can reproduce the results of conventional second order 



1 



perturbation theory to very high accuracy while using in- 
significant amounts of computer time. 

In future papers we will calculate higher order dia- 
grams while providing a development of renormalization 
theory that will allow a straightforward automation of the 
calculation of perturbative quantities. We believe that the 
Sine function representation will be useful in a wide range 
of applications from the evaluation of experimental results 
to the easy prototyping of model theories. 

2 Sine Functions 

An overview of the properties of Sine functions relevant to 
Feynman diagram calculations was given in [^j , and ^. [Tc| 
give more thorough discussions. Here we simply restate 
a few important results. A generalized Sine function is 
defined by 

sin [tv(x — kh)/h] 



S k (h, x) 



(1) 



7r(x — kh)/h 

Using an integral representation of this function, it is easy 
to show that any function f(z) which is analytical on a 
strip of a complex plane centered around the real axis and 
with width 2d, has the following approximation 

oo -i-oo 

p oo » oo 

J f(z)dz* J dzf(kh)S k (h,z) = h f(kh). 



The error of this approximation 
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satisfies 



|A/(/i)|<C 7 



fe=-c 



-Trd/h 



(2) 



(3) 



(4) 



2 sinh(7r<i//i) 

Here C is a number independent of h. This formula shows 
that if h < d, the difference between the exact and the 
approximate results decreases exponentially with h. 

When we calculate the Sine function representation 
there are two sources of error. Firstly, the quantity we 
are calculating differs from the exact value by an amount 
that depends on h, and secondly we are performing a nu- 
merical evaluation of infinite sum, and doing so inevitably 
introduces truncation errors and inaccuracies due to the 
finite precision of computer arithmetic. Consequently, by 
choosing h we attempt to establish the amount of accu- 
racy we hope to obtain from our calculation, but we must 
ensure that our numerical evaluation of the resulting sum 
is carried out to a commensurate level of precision. 



3 Propagators 
3.1 Coordinate Space 

In this section we obtain Sine function representations 
of the fermion and vector propagators in coordinate and 
momentum space.^j In order to achieve this, an approach 
similar to the one used in j(| will be applied to calcu- 
late coordinate space propagators. The momentum space 
propagators will be obtained by Fourier transformation. 

We use the fermion propagator to illustrate this cal- 
culation. Almost identical procedures are used to obtain 
expressions for the other propagators. 

In Euclidean space, the fermion propagator can be 
written as 



S F {x) = 



(2tt) 4 p 2 +m 2 ' 



(5) 



After introducing a cut-off and exponentiating the denom- 
inator we obtain: 



S F (x,A) = 



ds 1 



(to - 7> M ) x 



exp 



ipx 



(6) 



The second integral in ([|) is a Gaussian, which can be 
easily computed analytically: 



S F (x,A) 



ds e~ 



1 



2to ( s + 



32^ + 
1 



A 2 



exp 



(7) 



To make s dimensionless we shift s to s/M 2 , where 
M is an arbitrary constant with dimensions of mass. The 
final result does not depend on a specific value of M. 



2 In this paper we only directly consider massive, spin | fermions, 
but will often ignore this qualification when referring to the fermion 
propagator. However, these calculations will generalize to massless 
spin-^ fields and spin-3/2 fields. Likewise we will use "photon" 
interchangeably with "vector" since our explicit examples are drawn 
from QED, although we can compute diagrams containing massive 
gauge fields. 
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After performing this transformation and setting s to be 
equal to e 2 we obtain: 



S F (x,A) 




(8) 



At this point we use the Sine expansion (g) to arrive at the 
final form of Sine function approximation to the fermion 
propagator in coordinate space. 



S Fh (x,A) 
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In the above expression, and in the rest of this paper, we 
will use the subscript h to differentiate between the exact 
propagator and its Sine function approximation. We also 
introduce a shorthand notation 



a(ki,k 2 , k 3 , ...) 



= e klh 



,k 2 h 



?k 3 h 



M 2 , s 

(10) 



which makes algebraic expressions similar to (|^) more 
compact. 

By repeating all of the above steps, we find an approx- 
imation for the scalar field propagator: 



G h (x,A) 



M 2 h 
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exp 
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m 2 M 2 x 2 
M 2 ~ 4a(fci) 



(11) 



This form of scalar propagator becomes identical to the 
form derived in |fj), if M is set equal torn. In general, the 
numerical value of the parameter M does not have any 
significance for either scalar or fermion field propagators. 

Previously, we avoided introducing undetermined pa- 
rameters in the theory by shifting s to s/m 2 instead of 
s/M 2 . However, in the case of the photon, m = 0, so 
we cannot set M equal to the particle mass without di- 
viding by zero. Moreover, in order to regulate infrared 



divergences we will find it is necessary to assign a small 
mass, m v , to the photon and then take the limit where 
this mass goes to zero. In this case shifting s to s/m 2 
would produce terms involving m u /A, where m v — > and 
A — > 00. We prefer not do this, since this choice prevents 
us from handling ultraviolet and infrared divergences in- 
dependently. 

In an arbitrary gauge, the Euclidean space photon 
propagator is 

d%p 1 



D» v (x) 



(- 



(2tt) 4 p 2 + 
S^p 2 + (1 



1 

m 2 p 2 



(12) 



There are two factors containing p 2 in the denominator of 
(|l2|). The Sine expanded form of this expression requires 
independent exponentiation of each of these factors and 
consequently there are two independent summation pa- 
rameters, fci and k 2 , in the expression for the propagator: 
M 2 h 2 x e kih e k 2 h 

X 
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(4tt) 2 ^ 
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a(k 1: k 2 ) 
x 2 M 2 



k 2 h 



4a(fci, k 2 
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M 2 



x^x") 



(13) 



A.a{k\,k 2 ) 

The double summation in ( jl3| ) makes the task of eval- 
uating Feynman diagrams containing many photon lines 
significantly more difficult. Fortunately it is possible to 
solve this problem by choosing a particular gauge. In the 
Feynman gauge (£ = 1) the photon propagator is 



D^(x) 



(2tt) 4 p 2 



mf 



This can be approximated by 



D^(x,A) 



5" v 



+ OO 



M 2 h 



(4tt) 2 ^ 

v ; fel=-O0 
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(14) 



Hh)) 2 

M 2 x 2 



exp i- e ^M^~Mfco;- (15) 

At this point we have completed the derivation of Sine 
expanded versions of the propagators in coordinate space. 
However, for a number of reasons it is more convenient to 
evaluate diagrams in momentum space. The transition 
to momentum space significantly simplifies the algebraic 
form of the Sine function representations of the propaga- 
tors. 
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3.2 Momentum Space 

One way of obtaining approximate versions of the prop- 
agators in momentum space is to take Fourier transfor- 
mations of the corresponding expressions in coordinate 
space. Alternatively we can start by Sine expanding exact 
propagators in momentum space and then use the integral 
representation of the Sine function to get the final result. 
We use the former approach since it involves less algebra. 

The Sine function representation of the fermion field 
propagator can be written 



hM 2 ^ e k ^ h 

ki— — oo 



d A x [2ma(h) - iM 2 ^x^\ x 
m 2 x 2 M 2 



'-MM - ekl '\ T 2 



M 2 4a(fci) 



ipx (16) 



After taking the Gaussian integral and simplifying we 
obtain: 



S Fh (p, A) 



M 2 



+°° / 2 \ 

E ^fa-^w) 



ki— — oo 

( m - 7 ^)expf-^a(fci)). (17) 

Applying the same procedure to the scalar propagator 
(0) gives 



+ 00 



G h (p,A) 



M 2 



ki — — oo 



expik 1 h-e klh 



M 2 



exp[-j^a(h)). (18) 



Finally, the general photon propagator ( |l3| ) and Feynman 
gauge photon propagator (na) become 



h 2 exp((fc 1 +fc 2 )/ l -e fc ^^) 
M 2 ^ 



a(fci, fc 2 ) 

a(ki, k 2 ) 



(19) 



and 



S^JL y exp (k 1 h-ef tlh ^) 
M 2 I M 2 

fci=-oo 



exp 



(20) 



All the machinery necessary for us to start evaluating 
Feynman diagrams is now in place. The simplicity of the 
propagator in the Feynman gauge causes us to favor this 
form of the photon propagator (|2^) over the more gen- 
eral form (19). The disadvantage of this approach is that 
( p0| ) is not explicitly transverse, unlike when calculating 
in the Landau gauge, £ = 0, which automatically ensures 
that all computations are consistent with current conser- 
vation. Consequently, using the Feynman gauge requires 
us take extra steps to make sure that, despite approxima- 
tions, the current remains conserved. We will address this 
issue in detail when evaluating representative diagrams. 
However, before calculating the Feynman graphs we need 
to determine the error introduced by Sine expanding the 
propagators. 

3.3 Numerical Evaluation 

Since the momentum space propagators were obtained by 
manipulating coordinate space versions of the same prop- 
agators the error in the Sine function approximation to 
the momentum space propagator is very similar to that 
of the position space version. This can be verified by 
direct calculation. We work in momentum space in all 
further calculations. 

It is convenient to start by comparing the approximate 
version ( p7| ) of the fermion field propagator with the exact 
propagator 

m - 7 M p M 



Sf(p) 



P 



(21) 



Both ( |17| ) and ( J2 1| ) contain the matrix factor m — 7 M p M , 
which does not depend on the summation parameter k\ 
and thus cannot affect the accuracy of the approxima- 
tion. This allows us to define a function &(p,h), which 
estimates the error associated with the Sine expansion of 
the propagator. 



90, h) 



h 
M 2 

exp 



E 

fcl — — OO 



exp ( k\h — e 



ki h 



M 2 



-h a{ki) 



p 2 + m 2 



(22) 



It is easy to show that the same function can be applied 
to calculate the accuracy of the approximations of the 
scalar field propagator (|lj) and the photon field propaga- 
tor @. 

Figure |l| shows &(p, h) as a function of h. The results 
presented in this plot were obtained using 64 bit precision 
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48-14 r 




-1.2e-13 1 1 1 1 1 ' 1 

0.24 0.25 0.26 0.27 0.28 0.29 0.3 

Figure 1: Q(p,h) as a function of h is shown for m = 
0.511, M — 0.511, p = 0.75 and A — » oo. Computed using 
64 bit precision arithmetic. 



3B-06 r 




-5e-06 1 1 ' 1 1 1 ' ' I 1 

0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 



Figure 2: 0(p, ft,) as a function of ft is shown for m — 
0.511, M = 0.511, p = 0.75 and A -> oo. Computed using 
32 bit precision arithmetic. 



arithmetic. From the graph, it is obvious that accuracy 
better than 1 part in 10~ 14 can be achieved by taking the 
parameter ft, smaller than 0.28. It is also clear that the 
approximation stops improving after ft, gets smaller than 
0.25. At this point value of 0(p, ft) becomes smaller than 
the numerical noise produced due to the finite precision 
of the computation. 

To make sure that the above explanation is valid and 
that no other factors limit the accuracy of the approx- 
imation, 0(p, ft) was recomputed using 32 bit precision 
arithmetic. The results of this calculation are presented 
in Figure |^. The smallest error in this case is around 1 
part in 10~ 6 , which corresponds to ft, = 0.5. The behav- 
ior of G(p, ft) is dominated by numerical noise when ft is 
smaller than 0.5. 

The fact that the value of ft at which the approxima- 
tion breaks down increases by a factor of 2 as the arith- 
metic precision of the computation decreases by the same 
factor is a consequence of the the following result, derived 
in§: 

\Q(p,h)\<Ce-V h . (23) 

The error associated with calculations performed at finite 
precision is proportional to 10 _JV , where N represents the 
number of significant digits. By comparing these two re- 
lationships, we conclude that increasing N by some factor 
A should expand the range of ft in which calculations can 
be performed with adequate precision to include values 
up to \, 

The speculation presented above, can be verified 



directly. Trivial manipulation produces the following re- 
sult, 



1 



| log | &{p,h)\ 



< 



hb-d 



— as ft — » 0, 
a 



(24) 



where 6 is a constant that is related to C. 

Figure || shows that for sufficiently small values of ft all 
maxima of the function 1/| log 8(p, ft)|| fall on a straight 
line. This graph can also be used to determine the small- 
est value of the parameter h at which the accuracy of the 
approximation is not limited by the numerical precision 
of the calculation. 

Finally, it is necessary to make sure that the approxi- 
mation holds for any momentum. A plot of (p 2 +m 2 )Q(p, h) 
is shown in Figure |]. Here we see that by setting h to be 
equal to 0.275 it is possible to achieve a relative accuracy 
of the order of 1 part in 10~ 14 for a wide range of different 
momenta. Further, rapid oscillations of the error allow us 
to take integrals over momentum without a significant loss 
of precision. 

4 Evaluating Feynman Diagrams 

4.1 The Electron Vertex Function 

We use the Sine function representation to compute one- 
loop diagrams in QED, in accordance with the Sine func- 
tion Feynman rules, given in Since results for these 
diagrams are well known, we can use them as a first test of 
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Figure 3: l/ \ log |0(p, K)\\ as a function of h is shown for 
to = 0.511, M = 0.511, p = 0.2 and A -> oo . 



0.2 0.4 0.6 0.8 



1.2 1.4 1.6 1.E 



Figure 5: The Electron Vertex Diagram 



the approximation. Furthermore, it will be necessary to 
have expressions for these graphs to perform higher order 
calculations. 

We start by evaluating the electron vertex. This dia- 
gram, shown in Figure]^, contributes a/ (2tt) to the anoma- 
lous magnetic moment of the electron. In the following, 
we reproduce this result using the Sine function represen- 
tation. 3 

The expression for the leading correction to vertex 
function is 



u{p')5T»{p,p')u(p)=u{p') 



2tt) 4 

(-e lv )S F (k')^S F (k)(-e lp ) 



1 D pl, (p~k)x 

u(p). (25) 



After substituting our expressions for the propagators 
we find: 



u(p')5T»(p,p')u(p) = 



nip') 



2 h 3 



M 6 (2tt) 



E 

k1.kn.k3 



d%k x 



exp <j (fa +k 2 + k 3 )h - e klh m " 



2 

V 



k 2 h , k 3 h\ m e 



7p (m e - (k a + q a ) la ))Y l (m e - k x lx )Y x 



exp 



(P ~ kf 
M 2 



a{k\ 



(k + gf 
M 2 



■a(k 2 



Figure 4: Relative error of the approximation is shown 
by plotting (p 2 + m 2 )Q(p, h) as a function of p for to = 
0.511, M = 0.511, h = 0.275 and A -> 00. 



3 In this paper we focus our attention on numerical calculations 
of graphs which have been evaluated analytically, to ensure we have 
benchmarks for assessing the accuracy of the Sine function repre- 
sentation. The analytic results quoted in this section can he f ound 
in any thorough QED text, such as Peskin and Schroeder pi). 



G 



k 2 



*(p). 



(26) 



Table 1: Correction to the magnetic moment for different 
values of h. 



We eliminate gamma matrices with repeated indices, 
take the Gaussian integral over k, use the Dirac equation 
and 



l^P^uip) = -m e u(p), 
u(p')7V, t = -u(p)m e 



(27) 
(28) 



to simplify the result. Performing all the steps described 
above and using the Gordon identity, 



u(p')^{p) = u(p') 



p'» + p^ ia^qv 



2m„ 



2m. 



u(p), (29) 



to eliminate terms proportional (p' 1 +p' M ) we obtain the 
following expression for the second order correction to the 
electron magnetic moment: 
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Correction x — 


Error x — 
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0.4999999999949221 
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0.5000000000008583 
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0.499999900285 
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0.499998 


2.3 x 10~ B 
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1 a(ki 



7T M 2 ^ A 2 A 
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Figure 6: Electron self-energy diagram. 



„k 2 h 



M 2 



exp 



A/ 2 A 



(a(fc 2 ) + a(fc 3 ))a(fci) 



| (30) m = ~ J -^{-^)S{p-k){-e lv )D^{k) (31) 



Here A = a{k x ) + a(k 2 ) + a(k 3 ). 

The results of numerically evaluating this sum for dif- 
ferent values of h and m e = M = 0.511 are presented 
in Table |l|. As h is lowered the error decreases to a few 
parts in 10 13 . However, when h < 0.3 the accuracy de- 
creases. This is due to the use of fixed precision (64 bit) 
arithmetic, and is not intrinsic to the Sine function rep- 
resentation itself. 

Computationally, the triple sum is simple to evaluate, 
and for typical values of h it can be computed on a desktop 
computer in less than a second. In general, the time taken 
to evaluate a given diagram is a function of the number 
of internal propagators. The sums generated from QED 
diagrams are not significantly more complex than those 
derived from scalar diagrams with the same number of 
propagators. More detailed analysis of the relationships 
between the form of the diagram and the complexity of 
its Sine function representation, as well as discussions of 
efficient summation algorithms, are given in j| 

4.2 Electron Self-Energy 

The second order electron self-energy diagram is shown in 
Figure |[ The quantity, which represents a shift of 

the electron mass, can be written as follows. 



Substitution of the Sine function representations of the 
propagators transforms this to 



2„2 



h 2 e 



(2?r) 4 M 4 

exp j(fc 1 + fc 2 )h-e fel ' 1 



k-i .k 



M 2 



„k 2 h ■ 



d 4 k( ltI (m e - ^pc + ^k a )Y) x 



exp 



(P'k) 2 , k 2 



M 2 



M 2 



M 2 J 



(32) 



After taking the Gaussian integral, eliminating gamma 
matrices and simplifying, we obtain 



S(P) 



2m e +p° la (l-^±) 



fel >&2 



exp <j (h + k 2 )h - e fclh 4^ - e k - h ^ > x 



A 2 
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M 2 



exp 
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M 2 A 



(33) 



Here A = a(k%) + a(k 2 ). 
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0.3 0.31 0.32 0.33 0.34 0.35 0.36 0.37 0.38 

Figure 7: The absolute error in T,n(p) as a function of h 
is shown for m — 0.511, M — 0.511, p = 0.75 and A — > oo. 

The quantity above diverges as A goes to infinity. In 
order to renormalize it, we follow standard procedures. 
First we isolate the matrix dependence by the decompo- 
sition, 

£( P )=A(p 2 )+K 7CT iV), (34) 

where A(p 2 ) and B(p 2 ) are easily determined from equa- 
tion d33|). The finite quantity is formed by sub- 
traction of the divergent parts of S(p). 

E H (p) = (A(p 2 ) - A{0))+p° lrJ {B{p 2 ) - B(0)). (35) 

Combining expressions (|35[) and (|33|) we obtain: 



S«(p) 



/o\ ft* ^ 2m e +p CT 7 a (l " ^s- 
\n) 2 ^ 



A 2 

,2 



exp 



a(ki)a(k 2 )p 2 
M 2 A 



- 1 



(36) 



The difference between the values of £.r(p) obtained 
from (|3^) and the exact propagator correction to this or- 
der, 



a\ h 



(?) 



log 



(2m e + xp' 7 7 (T ) x 



(1 — a;)™ 2 + xm 2 



(1 — .T)m 2 + xm 2 + x(l — x)p 2 



.(37) 



Figure 8: Electron loop diagram. 

is shown as a function of h in Figure |^. By setting h to 
be equal to 0.31 it is possible to achieve accuracy better 
than 1 part in 10 in this calculation. 

4.3 Charge Renormalization 

To second order in perturbation theory, the electron loop 
diagram IP 1 ' (p) , shown in Figure ||, is the only relevant 
graph. II A " / (p) has the familiar form: 



W{p) 



d 4 k 



tr [^S(k)YS(p + k)} . (38) 



After substituting Sine expanded representations for S(k) 
and S(p + k) and taking the trace, we obtain 



2 h 2 



e 2 h 



4Af 4 7T 4 
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fcl,fc 2 



exp < (ki + k 2 )h - (i 



M 2 



y d 4 k[ P ^k" + 2/c^ + fcv - 

• fc + A: 2 + m 2 )] x 
f k 2 (p + k) 2 
° XP | "M2 a( 2 > M*— a ( kl > 



(39) 



By now it is not surprising that the integration over 
internal loop momentum k can be easily performed ana- 
lytically. As we mentioned earlier, in order to use this (or 
any approximation) we must be sure that current conser- 
vation is respected. We achieve this by multiplying IP" 
by the transverse projection operator 



pXpa 
p2 



rArr 



(40) 



After taking the integral in (B9h and projecting out the 



transverse part with (40) we obtain 



n*-» = V£iex P { 



a{k 1 )a(k 2 )p 2 
AM 2 



8 



exp i (fcj + k 2 )h - (e felh + e fe2/i ) ^ > x 



a(fci)a(/c 2 )p 2 



It is convenient to make the usual definition: 



ip» = n( P 2 ) 



(42) 



The quantity II (p 2 ) diverges as A — ► oo. We can renor- 
malize it by subtracting out divergent parts, which are 
contained in the first two terms of the Taylor series 



n^ 2 ) = iiV)-n(o)-p 2 du{p2) 



dp 2 



(43) 



p 2 =0 



Combining the definitions above with equation (ffl|), we 
derive the final result 



7T — A 2 

fei ,k 2 



<>K V { (k 1+ k 2 )h~ ( e *l" +e tafc)^|J X 



a(fci)a(fc 2 )p 2 



exp 



A 2 

a(/ci)o;(A:2)p 2 



AM 2 



M 2 a(ki)a(k 2 ) 2 



P I A 



M 2 



(44) 



As before we check the accuracy of the above equation 
by comparing it to the exact result, 



dx x(l — x)p 2 x 



log 



m 2 + x(l — x)p 2 



(45) 



The result of this comparison is presented in Figure |^. 
The error of the approximation stops decreasing after it 
reaches a value close to a 1 part in 10 15 at h = 0.3. 

5 Conclusions 

In this paper we have extended the Sine function tech- 
niques developed in [|| to include fermion and vector prop- 
agators. We display approximate forms of scalar, fermion 
and photon field propagators derived in both coordinate 




0.31 0.32 



0.34 0.35 



Figure 9: Absolute error in Hji(p 2 ) is shown as a function 
of ft for m e = 0.511, M = 0.511, p= 0.75. 



and momentum space. The ease of obtaining high accu- 
racy with these approximations is demonstrated by per- 
forming several numerical comparisons between exact and 
Sine expanded versions of the propagators. These com- 
parisons show that using 64-bit numerical precision, cur- 
rently available on all hardware, it is possible to achieve 
errors as low as a few parts in 10 15 . On the other hand, if 
high precision can be sacrificed for speed, it is possible to 
increase the parameter h with a corresponding decrease 
in accuracy. However, even the highest accuracy calcu- 
lations described in this paper typically take less than a 
second on standard desktop computers. 

To illustrate the use of the Sine function representa- 
tion and lay the ground for future work we evaluated three 
one- loop QED Feynman diagrams. In all three cases it is 
easy to perform the calculations by hand, so we were able 
to evaluate the error in computations based on Sine ap- 
proximated propagators. In general, there are three fac- 
tors that limit the accuracy of the calculation. First, some 
error is introduced when we use the Sine representation 
to approximate the propagators. It is possible to improve 
the approximation by decreasing the value of h. Second, 
truncation of the sums, which have to be evaluated in or- 
der to compute Feynman graphs using the Sine expansion, 
produces another contribution to the total error. Finally, 
the finite numerical precision of the computation is the 
third source of error. 

The three sources of error are closely related. For ex- 
ample, an attempt to achieve higher accuracy by decreas- 
ing the value of h without simultaneously increasing the 
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range of the summation would result in a loss of preci- 
sion. To minimize the error, all parameters must be set 
to optimal values. While it is easy to manipulate h and 
the range of the summation, increasing numerical preci- 
sion can be difficult since we often have to evaluate terms 
which involve subtracting two almost equal quantities. 

The results in this paper provide further confirma- 
tion that the Sine function representation can be used 
in high precision calculations with the expenditure of sur- 
prisingly small amounts of computing time. Higher or- 
der calculations will require significantly more time than 
those presented in this paper, but results to date indi- 
cate that Sine techniques can still be orders of magnitude 
faster than Monte Carlo methods when high accuracy is 
required. Finally, the methods we have presented here 
can, in principle, be automated. In future work we plan 
to focus on this task in order to facilitate the calculation 
of higher order perturbative graphs in arbitrary theories. 
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